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SUMMARY 

A steady-state computer simulation model of an Electrochemical Depolarized 
Carbon Dioxide Concentrator (EDC) has been developed. The mathematical 
model combines EDC heat and mass balance equations with empirical corre- 
lations derived from experimental data to describe EDC performance as a 
function of the operating parameters involved. The model is capable of 
accurately predicting performance over EDC operating ranges. Model simula- 
tion results agree with the experimental data obtained over the prediction 
range . 

Eight experiments were conducted in support of the mathematical model to 
study the effect of the cell matrix on cell internal resistance (formation 
factor); anode gas pC0 2 levels; process air flow rate; current density, 
process air inlet pC0 2 ; matrix thickness; electrode thickness; thermal con- 
ditions within a cell; and the transient response of a cell to changes in 
process air inlet pC0 2 , current density and air flow rate. 

A computer program was written for the numerical solution of the mathematical 
expressions describing system performance. The program was written in 
Fortran IV for use on a Univac 1108 Digital Computer. 

INTRODUCTION 

In an Electrochemical Depolarized Carbon Dioxide Concentrator (EDC) numerous 
physical, chemical and electrochemical processes interact simultaneously. 

These processes can be experimentally defined and characterized mathematically 
for computerized solution to predict EDC performance as a function of the 
operating parameters involved. In addition to predicting performance, the 
mathematical model can be used to optimize subsystem designs and performance, 
and aid in integrating an Electrochemical Carbon Dioxide Concentrator Sub- 
system with other regenerative life support subsystems. 

In parallel to the development of a Six-Man Electrochemical Carbon Dioxide 
Concentrator Subsystem^ , a program was completed to 

1. define the electrochemical carbon dioxide (CO^) concentration 
process to analytically determine the relationship between EDC 
performance and the operating parameters and to provide a basis 
for assumptions made; 

2. perform eight experiments to supplement previously obtained 
data and to support the development of an analytical mathe- 
matical model, * 

3. transcribe the experimentally defined functional relationships 
between operating parameters and EDC performance into mathe- 
matical expressions (empirical correlations) for computerized 
solution; and 


(1) All references cited are listed at the end of the report. 


1 



JCife Systems. Juc. 


4. prepare a computer program to predict EDC performance by solving 
the empirical correlations obtained and EDC heat and mass bal- 
ance equations as a function of the operating parameters 
involved. 

The theoretical work completed was used to gain an understanding of the processes 
involved in the electrochemical concentration of C0 2 and to define where experi- 
mental results were needed. The computer program prepared was based on experi- 
mental performance results and did not necessarily incorporate the theoretical 
work completed. 

The following sections summarize the results of the program. 

ELECTROCHEMICAL CARBON DIOXIDE CONCENTRATION 

Carbon dioxide is removed from a flowing air stream as it passes over the cathode 
of an EDC cell. Each cell consists of two porous electrodes separated by a 
porous matrix containing an aqueous solution of cesium carbonate (Cs 2 C0_). 

Plates adjacent to the electrodes provide passageways for distribution of the 
process gases and electrical current over the electrode surfaces. The specific 
electrochemical and chemical reactions are detailed in Figure 1. 

Moist air containing C0 2 is fed into the cathode where the electrochemical 
reaction_of oxygen (0 2 ) in the air, water (H 2 0) and electrons forms hydroxyl 
ions (OH") . Hydroxyl ions then react with the C0 2 forming carbonate ions (CO ") . 
The output from the cathode compartment is moist air at a reduced CCL partial 3 
pressure (pC0 2 ), At the anode side hydrogen (H 2 ) is fed into the cell. The 
electrochemical reaction of the H 2 and OH to form H 2 0 and electrons decreases 
the concentration of OH in the electrolyte at the anode. Carbon dioxide is 
produced in the resulting equilibrium shift completing its transfer from the 
cathode compartment to the anode compartment. The output from the anode compart- 
ment is C0 2 mixed with unreacted The overall reaction is exothermic and 

is accompanied by the formation of electrical energy. 

EDC Performance 

The performance of an EDC is reflected by CCL removal efficiency and electrical 
efficiency. ^ 

Carbon Dioxide Removal Efficiency 

Inspection of the overall reaction as based on the C0„“ transfer mechanism 

0 2 + 2CC> 2 + 2H 2 *»■ 2C0 2 + 2H 2 0 + Electrical Energy + Heat 

shows that two moles of C0 2 can be transferred for one mole of 0_ consumed. 

This represents a C0 2 removal efficiency of 100%, The equivalent weight ratio 
is 2.75 lb of C0 2 removed for each pound of 0 2 .consumed. This ratio has been 
referred to as the Transfer Index (TI) . 
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START/STOP SWITCH 



CATHODE 

ELECTRODE 


ANODE 

ELECTRODE 


CATHODE REACTIONS : 

0 2 + 2H 2 0 + 4e" = 40H‘ 

40H~ + 2C0 2 = 2H 2 0 + 2C0 3 ' 


ANODE REACTIONS : 

2H ? + 40H~ = 4H 7 0 + 4e“ 

2C0 3 = + 2H 7 0 = 40H" + 2C0 2 


OVERALL REACTION: 


0 2 + 2C0 2 + 2H 2 


2CO ? + 2H ? 0 + Electrical Energy + Heat 


FIGURE 1 ELECTROCHEMICAL AND CHEMICAL REACTIONS 
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To gain a better understanding of the dependency of the C0 2 removal process 
and, hence, removal efficiency upon the operating parameters involved, the CO 
mass transfer process was subdivided into five steps. They are (1) air cavity 
transfer, (2) cathode reactions, (3) electrolyte transfer, (4) anode reactions, 
and (5) cavity transfer. 

Air Cavity Transfer . Carbon dioxide in the process air stream diffuses in 
the cathode compartment to the electrode surface and then through the elec- 
trode pores to the gas-electrolyte interface. The air cavity CO transfer 
rate can be estimated using a convective mass transfer coefficient. The 
equation for air cavity C0 2 transfer is 

N C0 = (k c )(A)(M CO ^ Cp C0 2 - PC0 2 °)/RT 
2 2 

Where 

N C0 = ^2 trans ^ errec ^> lb/hr 

k c - convective mass transfer coefficient, ft/hr 

2 

A = transfer area, in 

pCO^ = C0 2 partial pressure in air stream, psia 
pC0 2 ° = C0 2 partial pressure at the gas-electrolyte interface, psia 
M C0 = mo ^ ecular weight of C0 2 , lb/ lb mol 
R = Universal Gas Constant, ft- lb/lb mol-R 
T = absolute temperature, R 

Air cavity C0 2 mass transfer rate, therefore, depends on pCO • k , a function 
of process air flow rate: cathode compartment and electrode geometry and the 
physical properties and diffusion coefficients of the gases involved; tem- 
perature; and pressure. The physical properties and diffusion coefficients 
can be assumed constant with negligible error over typical EDC operating 
ranges in temperatures (60-80F) and pressure (14.7 psia). For a given EDC 
cell geometry and normal temperature and pressure conditions, air cavity CO 
mass transfer rate is assumed to be only a function of process air flow rate 
and pC0 2 . 

Cathode Reactions . Carbon dioxide reacts with electrochemically-generated 
OH at the cathode to form CO^ in two consecutive reactions^ ^ 

C0 2 + OH" HC0 3 “ 

hco 3 " + OH" co 3 = + h 2 o 

In strong basic solutions the equation for the reaction rate is^ 

-d [CO ] 

= k [C0 2 ] [OH - ] 


4 



Cifc Systems, Jhc. 


Where 

k = second order rate constant, 1/mol-sec 
[C0 o ] = CO^ concentration, mol/1 
[0H~] = OH - concentration, mol/1 

The C0 ? available for reaction is in equilibrium with pCO^ 0 at the gas- 
electrolyte interface in the electrode pores. The OH available for reaction 
is a function of cell current density, and k is a function of temperature. For 
the cathode reactions, CO^ transfer rate is a function of pCO^, cell current 
density and temperature. 

Electrolyte Transfer . Carbon dioxide as CO^” and unreacted OH diffuse 
through the bulk electrolyte to the anode. Since current is carried by CO^ , 

C0 9 mass transfer is the sum of molecular and convective diffusion and ionic 
migration under the electrical potential gradient between the two cell electrodes. 
The equation describing the transfer of charged species through electrolyte 
solutions is given by the Nernst-Planck equation* - } 


N. 

1 

Where 

D. = diffusion coefficient of j 
J 

Zj = electrical charge of j 

F = Faraday’s constant 

Uj = electrical mobility 

C. = concentration of j 
J 

0 = electrical potential 
x = transfer distance 
v = bulk fluid velocity 

For transfer through the electrolyte, v can be assumed to be zero. The mass 
transfer rate through the electrolyte then is a function of the diffusion co- 
efficients, electrolyte composition, matrix structure, temperature and current 
density. 

Anode Reactions. Carbonate ions react with electrochemically formed H^O to 
form CO As H 2 0 is formed the equilibrium shifts from CO," to bicarbonate 
ions (HCO^ ) . 

H 9 0 + C0 3 = HC0 3 - + OH" 


„ dC i „ „ d0 

= - D. j— z . Fu . C . + C . v 

J dx j 3 J dx j 
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At the lower anolyte pH levels, two parallel reaction mechanisms can take 
place 


1 . 

HC0 3 

~ — *~C0„ + OH 

2. 

hco 3 - 

+ H O ^T^TH CO 


JL h 2 o ♦ C0 2 

The equation describing the reaction rate is 
d[CO ] k 

— af— = [HC0 3 “] - kj [C0 2 ] [OH'3 + k 2 [H 2 C0 3 ] 

Where 

k^, k 2 * reaction rate constants 
K^, K 2 = equilibrium constants 

a H 2 0 = activity of H 2 0 

The concentration of dissolved CO 2 is in equilibrium with the anode gas pC0_ 
at the gas-electrolyte interface m the electrode pores. For the anode reactions, 
C0 2 transfer rate is a function of anode gas pC0 2 , cell current density and 
temperature. 

Hydrogen Cavity Transfer . Carbon dioxide diffuses through the electrode pores 
and into the anode gas stream. Hydrogen cavity C0 2 mass transfer rate can be 
estimated using a convective mass transfer coefficient similar to the air cavity 
mass transfer coefficient. The equation for H 2 cavity C0 2 transfer is 

N C0 = (k c )(A:i(M C0 )(P C0 2° - P C( V /RT 
2 2 

Where 


k 2 ^ cc y ^ a H 2 o^ 


pCO 2 , pC0 2 ° and k c are evaluated for the H 2 cavity conditions. 

Electrical Efficiency 

The electrical energy produced by the electrochemical reaction in the EDC is 
a function of the current density and the cell voltage. The theoretical open- 
circuit voltage is 1.23 volts. In practical applications, and with current 
flowing, cell voltages of less than 1.23 volts result. Electrical efficiency 
is therefore reflected by cell voltage with high cell voltage representing 
high electrical efficiency. 

Cell voltage is the sum of the reversible electrochemical reaction voltages 
(E r ) overvoltages (E Q ) and concentration polarizations (E c ) at each electrode 
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and the voltage loss due to the internal resistance of the electrolyte-filled 
matrix (E TR ) * The expressions for E„, E„ and E_ can be derived from thermo- 
dynamic considerations. However, E_ R depends on electrolyte matrix conductivity 
for a given matrix material and configuration (i.e., matrix compression, 
porosity, etc.). 

A general equation can be used tp. describe cell voltage as a function of cell 
current density and temperature 1 

i -i 

E = E D - (a Ini + bi + cln (— r — ) + d) 

R X L 

Where 

а, b, c and d = constants 

E = cell voltage, volt 
i = cell current density, ASF 
T = cell temperature, F 
i T = limiting current density, ASF 

The constants in the equation depend on temperature and the electrode-electro- 
lyte matrix-electrode structure and properties and can be found using cell 
performance data for a specific EDC design. The limiting current density 
depends on anode gas flow rate and pH^, and air flow rate and pC^* 

Operating Parameters Affecting EDC Performance 

Six major operating parameters affect EDC performance. They are 

1. pCC^ in the air supplied to the cell, 

2. current density at which the cell operates, 

3. air flow rate through the cathode compartment, 

4. pC0 9 in the H ? stream, 

5. H 0 flow rate through the anode compartment, and 

б. temperature level at which the cell operates. 

The effect of each parameter on CO2 removal and electrical efficiencies has 
been characterized experimentally. 

EDC Hardware 

The mathematical model developed is based on Life Systems* EDC cell design. 
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Figure 2 is a schematic of the concentrator cell. The cell consists of a 
cathode and an anode gas compartment, two current collectors with gas compart- 
ment spacers, two electrodes, and an electrolyte matrix. Figure 3 is a photo- 
graph of the cell parts. The basic materials of construction are Life 
Systems’ baseline electrodes, asbestos electrolyte matrix, nickel expanded 
metal (Exmet) gas compartment spacers, nickel current collectors, and poly- 
sulfone structural cell frames. The cell utilizes externally air-cooled fins 
attached to the anode current collector to remove the waste heat generated 
by the electrochemical reactions. The materials of construction and cell 
dimensions are summarized in Table 1. 

EDC Control 

Steady-state EDC operation requires that H-0 and heat generated in the electro- 
chemical cell reactions be removed to maintain cell moisture balance. Cell 
temperature is controlled at a set differential temperature above the process 
air inlet dew point temperature. This control scheme allows the process air 
humidity and cell electrolyte concentration to remain in balance as generated 
H^O and heat are removed. 

Water Removal 


The H^O generated in the electrochemical reactions must be removed from the 
cell by the process air stream. (Some H^O is removed in the stream but 
due to the low flow rates this can be assumed negligible.) For steady-state 
it is necessary to balance the amount of HO that is added to the process air 
stream such that the electrolyte concentration remains constant. The equation 
for the H^O produced in an EDC is 

W = 7.4 x 10" 4 (N) (I) 


Where 


W = H_0 produced, Ib/hr 
N = number of series-connected cells 
I = cell current, amp 

Process air humidity conditions in which all the H^O is not removed would 
result in a decrease in the electrolyte concentration. The resulting volume 
increase may cause electrolyte to be deposited in the gas cavities. Should 
more than the electrochemical ly-formed H^O be removed, the electrolyte con- 
centration would increase and the resulting volume decrease may cause electrolyte 
precipitation or electrolyte to be pulled out of the electrode pores (the result 
is a cell failure allowing gas to cross from one cavity to another) . 

Heat Removal 


For the electrochemical reactions involved, the cell voltage derived from the 
heat of reaction, is 1.25 volt. The net heat produced by an EDC is given by the 
following equation. 
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FIGURE 2 EDC CELL SCHEMATIC 
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FIGURE 3 EDC CELL PARTS 
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TABLE 1 EDC CELL CHARACTERISTICS 


Active Area (per cell) 
Air Cavity Height 
Cavity Height 
Cell Matrix 
Cell Frames 
H 2 Current Collector 
0 2 Current Collector 
Gaskets 


0.244 Ft 2 (4.1 x 8.6 In) 
0.082 In 
0.062 In 

0.030 In, Asbestos 
Injection-Molded Polysulfone 
Nickel 200 
Nickel 200 
Ethylene Propylene 
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Q = (1.25-E) NI x 3.419 

Where 

Q = heat produced, Btu/hr 
I = the current, amp 
E = cell voltage, volt 

The net heat produced is absorbed by the process air, the H. stream and the 
cooling air so that 1 

Q = AH a + + AH c 

Where 

AH^ = heat removed by the process air 

AH^ = heat removed by the process H 2 

AH^ = heat removed by the process cooling air 

TEST PROGRAM 

A successful Test Program supplemented previously obtained performance data for 
1-, 3-, 15- and 90-cell EDC designs and produced specific operating data for 
current density, pCC> 2 and air flow variations that would aid in mathematically 
characterizing EDC performance. 

The testing was parametric in nature and a single-cell EDC was used. The 1 
three-cell test facility also developed under NASA Contract NAS2-6478 was ' 
modified for single-cell testing and was used for the experiments. A schematic 
of the test facility is presented in Figure 4. Specifically, the test 
facility has the capability of varying process air flow rate, dew point and 
pC0 2 , H 2 flow rate, current density and cell temperature. 

A series of eight experiments were completed to support the EDC mathematical 
model development. They were 

1. formation factor determination, 

2. effect of anode gas pC0 2 variations, 

3. air cavity mass transfer coefficient determination, 

4. effect of low current density, 

5. effect of matrix thickness, 

6. effect of electrode thickness, 
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7. EDC temperature profile determination (thermal gradient), and 

8. EDC transient response to step changes in operating conditions. 

The baseline test conditions of the experiment are given in Table 2. Table 
3 summarizes the ranges over which the test parameters were varied. 

Formation Factor Determination 

The formation factor is defined as the ratio of the observed internal resistance 
(IR) of the cell (electrolyte-filled matrix) divided by the IR of the electrolyte. 
The formation factor for Life Systems’ EDC cell design was determined for two 
concentrations of Cs-CO,. A Kordesch-Marko (K-M) Bridge was electrically con- 
nected in series witn a single-cell EDC. The interrupted current from the 
K-M Bridge allowed the measurement of cell IR as a function of current density. 
The formation factor was calculated using the K-M Bridge IR measurements and 
Cs^CO^ conductivity data. 

The results of the current density span for each Cs^CO^ concentration are 
presented in Figure 5. The formation factor for the 61.5% Cs^CO^ charge con- 
centration was 6.2. For the 55% Cs^O^ charge concentration the' 5 format ion 
factor was found to be 5.8. (The formation factor calculations are given in 
Appendix A.) The results of the two formation factor determinations agree with 
experimental error and the average value, 6.0, can be used for the formation 
factor. 


Effect of Anode Gas pC0 2 Variations 

The effect of anode gas pC0 2 levels on EDC CO^ removal and electrical efficiency 
was determined for a pCO^ range of 0-690 mm Hg (at a total anode gas pressure 
of 760 mm Hg) . Carbon dioxide was added to the inlet stream to simulate 
the various pC0 2 levels. The results of the experiment are presented in Figure 
6. Carbon dioxide removal efficiency remained constant (TI = 1.98) through 
640 mm Hg while cell voltage decreased steadily at approximately 0.2 millivolts 
per mm Hg. At 690 ram Hg or 90.8% by volume of C0 2 in the anode compartment, both 
TI and cell voltage decreased to 0 and 0.09 volts, respectively. 

The minimum H 2 flow rate for normal EDC operation is 1.3 times the stoichio- 
metric amount required by current. The maximum EDC anode gas pCCL level then 
is 585 mm Hg (assuming 1.3 x stoichiometric H 2 flow and TI = 2 . 75 j . For normal 
operating conditions the anode gas pC0 2 level is never greater than 253 mm 
Hg (assuming 3 x stoichiometric H„ flow and TI - 2.75). Anode gas pC0„ level 
then has a negligible effect on EDC performance causing only a slight decrease 
in cell voltage over the EDC's typical operating range. 

Air Cavity Mass Transfer Coefficient Determination 

The air cavity convective mass transfer coefficient was determined for an EDC. 

The convective mass transfer coefficient, k , was estimated using the following 
equation (see Appendix B for derivation) 

k = AV b 
c 
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TABLE 2 BASELINE CONDITIONS FOR MATH 


Process Air Inlet 

Flow Rate 
Pressure 
Dew Point 

pco 9 

po 2 ‘ 

Process Inlet 

Flow Rate 
Pressure 


Cell 


Current (Current Density) 
Temperature 

Electrolyte 
Charge Concentration 


MODEL EXPERIMENTS 


0,44-0.50 Scfm 
Ambient 
50 ±2F 

3.0 ±0.2 mm Hg 

3.1 Psia 


300 ±30 cc/Min 
Ambient 


4.88 Amp (20 Asf) 
(Process air inlet 
dew point + 20) ±1F 
Cs ? C0 

6175 50.5% w/w 
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TABLE 3 SUMMARY OF MATH MODEL PARAMETRIC CHARACTERIZATION CURVES 


Experiment 

No. 

Description 

Performance 
Variable (s) 

Range ^ 

Curve (s) Reflecting 
Results 

1 

Matrix Formation 
Factor 

Current Density; 
Charge Concen. 

10-40 Asf for 55% and 
61.5% 

Figure 5 

2 

Anode Gas pC0 2 
Variations 

Inlet anode gas 
P C0 2 


Figure 6 

3 

Air Cavity Trans- 
fer Coefficient 

Air Flow; pCO^ 

0.2-0. 8 Scfm for 1.0, 1.5, 
2.0 and 10.2 mm Hg 

Figure 7 

4 

Low Current 
Density Test 

Current Density, 

pco 2 

5-30 Asf for 0,5, 3.0 and 
9.0 mm Hg (nominal) 

Figures 8 and 9 

5 

Thin Matrix Test 

Current Density 

10-40 Asf 

Figure 10 

6 

Electrode Thickness 

Current Density 

10-40 Asf 

Figure 11 

7 

Thermal Gradient 

Current Density, 
Inlet Relative 
Humidity 

80% and 63% RH at 20 Asf, 
80% RH at 40 Asf 

Figures 13 through 15 

8 

Transient Response 

Operating Time 
after step changes 
in Current, Density, 
Air Flow and pCO^ 

0-24 Hr for 20-30 Asf. 

30-20 Asf, 0.48-0.91 Scfm. 
0.91-0.49 Scfm, 3. 0-5. 5 mm Hg, 
5. 0-2. 8 mm Hg 

Figures 16 through 21 

Xa) Variable ranges not specified 

are at baseline conditions. 
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Cell Voltage. Volt 
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Where 


A, b = constants 

V = volumetric air flow rate 

At low inlet pC0 2 levels (-1.5 mm Hg at 20 ASF) air cavity mass transfer is 
rate-controlling and C0_ removal is a function of the process air flow rate 
and pC0 2 . At high pCO ^levels (-5.0 mm Hg at 20 ASF) air cavity mass transfer 
is not rate-controlling and, hence, air flow rate has less effect on CO^ removal* 
By measuring TI as a function of air flow rate at pC0 2 levels of 1.0, 1.5 and 
2 . 0 mm Hg, A and b were calculated giving k as a function of process air 
flow rate. Also, the air flow test was run for 10.2 mm Hg to show that air 
flow had little effect on TI at higher pC0 2 levels. The results of the experi- 
ment are shown in Figure 7. The calculated values of A and b for 1.0, 1.5 
and 2.0 mm Hg are 

A = 113.9, b = 0.49 at 1.0 mm Hg 

A = 90.8, b = 0.36 at 1.5 mm Hg 

A = 88.3, b = 0.41 at 2.0 mm Hg 

The correlation of data to the general equation derived for k c was good, 
verifying that air cavity mass transport is rate-controlling at -1.5 mm Hg. 

Effect of Low Current Density 

At low current densities (-10 ASF) only a small amount of CO^ is transferred 
making normal experimental transfer measurements inaccurate. For this experi- 
ment, a different technique was used to measure CO^ transfer. A tube of 
ascarite was placed in the anode gas exit stream such that the H^ and CO^ 
flow entering the ascarite tube and the H^ flow leaving the ascarite tube 
could be measured. Since all the C0 2 in the anode gas stream was absorbed by 
the ascarite, the difference in the flow rates of the two streams equalled the 
amount of CC^ transferred. 

EDC performance was determined as a function of current density for 5.0, 7.5, 
10.0, 20.0 and 30.0 ASF at 0.5, 3.0 and 9.0 mm Hg (nominal) inlet pCCL levels. 

The results of the experiment are presented in Figure 8. These results were 
cross-plotted in Figure 9 to show the effect of current density on EDC per- 
formance at a given pCCL level. The results of the experiment for 20 and 30 
ASF agreed with results previously obtained for TI as a function of current 
density and pCO^. 


Effect of Matrix Thickness 

EDC performance was determined as a function of current density for a 0.010- 
inch thick asbestos matrix as opposed to the standard 0.030-inch matrix used 
in Life Systems’ EDC cell design. The results of the experiment are presented 
in Figure 10. The C0 2 removal efficiency was lower for the thinner matrix 
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FIGURE 7 EFFECT OF AIR FLOW RATE AT FOUR PC0 2 LEVELS 
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FIGURE 8 EFFECT OF LOW CU 
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(TI = 1.2 as opposed to 2.1 for the 0. 030-inch baseline matrix) indicating that 
for the shorter diffusion length, cell concentration gradients were such that 
C0 2 backdiffusion inhibited CC> 2 transfer. Cell voltage was not affected. 

Effect of Electrode Thickness 

EDC performance was determined as a function of current density for 0.005 inch 
thick electrodes, as opposed to the standard 0.010 inch thick electrodes used 
in Life Systems’ EDC cell design. The results of the current density span are 
presented in Figure 11. The thinner electrodes increased TI from 2.1 for base- 
line configuration to 2.35. Cell voltage was unaffected. The increase in TI 
indicates that the electrode thickness improved air cavity convective mass 
transfer by decreasing the distance C0 2 has to diffuse through the electrode 
pores. 


EDC Temperature Profile Determination (Thermal Gradient) 

A temperature profile for an EDC cell was obtained for three process air heat 
load conditions: an air inlet relative humidity (RH) of 80% and 63% at 20 
ASF and 80% at 40 ASF. Ten thermocouples were placed in the air cavity adjacent 
to, but not touching, the cathode. Two thermocouples were inserted, one each 
in the air inlet and air outlet manifolds. Two thermocouples were inserted 
in the cell matrix in the same plane as two air cavity thermocouples. The 
thermocouple location schematic is presented in Figure 12. The temperature 
profiles for the three process air heat load conditions are presented in Figures 
13, 14 and 15. There was only a small difference in temperature between the 
thermocouples located in the cell matrix and those located in the air cavity. 

The air cavity temperature then reflects cell matrix temperature and can be 
used for moisture balance control. 

Transient Response 

EDC performance was measured as a function of time for step changes in current 
density, air flow rate and pC0 2 . Step changes were made for one increase and 
one decrease in each of the three operating parameters and EDC performance was 
monitored as a function of time until equilibrium was reached. The transient 
responses in EDC performance for the six step changes are presented- in Figures 
16 through 21. 

Transfer Index was measured on both the side and air side to determine 
transient differences in the C0 2 removed from the process air and concentrated 
into the process H 2 . No transient differences were observed as the EL side 
and air side TI measurements agreed in experimental accuracy. 1 

The step decrease in current density required the longest time for the EDC to 
regain equilibrium - 10 to 24 hours. The step increase in air flow required 
3 to 4 hours for the EDC to regain equilibrium while both pC0 2 changes and the 
increase in current density and decrease in air flow required approximately 
2 hours. Each cell voltage responded in 1 to 2 hours after a change in current 
density with only negligible effects being observed for the other parameter 
step changes. 
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FIGURE 15 TEMPERATURE PROFILE FOR 80% AIR INLET RELATIVE HUMIDITY, 40 ASF 
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FIGURE 16 TRANSIENT RESPONSE OF A CO CONCENTRATOR TO A STEP INCREASE IN CURRENT DENSITY (20 ASF TO 30 ASF) 
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FIGURE 17 TRANSIENT RESPONSE OF A CO CONCENTRATOR TO A STEP DECREASE IN CURRENT DENSITY (30 ASF TO 20 ASF) 
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FIGURE 19 TRANSIENT RESPONSE OF A C0 2 CONCENTRATOR TO A STEP DECREASE IN AIR FLOW CO. 91 SCFM TO 0.49 SCFM) 
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EMPIRICAL CORRELATIONS 

Experimental results graphically define the interaction between the major 
operating parameters and EDC performance. These graphically defined relation- 
ships were transformed into empirical correlations or mathematical expressions 
for use in the computer simulation program. 

Empirical correlations were required for the computer simulation model in four 
areas 

1. CO^ removal efficiency, 

2. cell voltage, 

3. moisture tolerance, and 

4. heat removal. 

CO^ Removal Efficiency 

The TI was experimentally found to be a function of process air inlet pCCL, 
current density and process air flow rate as predicted by theory. A combination 
of the results obtained in the mathematical model experiments and data selected 
from previous tests on 1-, 3-, 15-, and 90-cell EDC designs yielded the set 
of data points shown in Figure 22. These data provided ithe basis for the 
empirical correlation for TI as a function of process air inlet pC0~ and current 
density. A fifth-degree polynomial was fitted by the least-squaresmethod 
through the data points obtained for 10, 20, 30 and 40 ASF. The equation used 
in the data fitting program was 
5 

TI = L S. (pCO ) X 
i = 1 1 1 

Where 

S i = the constants for a given current density. 

Similar equations were added for 15, 25 and 35 ASF to aid in interpolating 
between curves. The values for S. are given in Table 4. Simultation curves - 
for 10, 20, 30 and 40 ASF drawn through the data points on Figure 22 verify 
prediction. As can be seen, the simulation curves obtained using the polynomial 
relationship describes the data within experimental error. 

The effect of process air flow was incorporated into the TI correlation by 
using a normalized value for the inlet pCCL as a function of air flow rate. 

The correlation derived for the normalized pCO_ was based on experimental results 
and is given by the following equation 1 

P = pCO x (V/0.44) ^°* 4S " °- 0045 pC0 2 2 ) 

& z 

Where 
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FIGURE 22 TI AS A FUNCTION OF pC0 2 AND CURRENT DENSITY 
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TABLE 4 VALUES FOR CONSTANTS 


Current Density 



Asf 

Si 


s 2 

S3 

S 4 




10 

1.75187 


-5 . 11449xl0 -1 

7.75073X10" 2 

-5 ,92584xl0" 3 

1.79862xl0" 4 


15 

1.42306 


-2.71935X10" 1 

2.10187xl0" 2 

-3.4839xl0" 4 

-1. 8475xl0~ 5 


20 

1.08190 


-1 . 1082 lxlO" 1 

-8.1408xl0' 3 

1.96359x10"^ 

-8.5832xl0" 5 


25 

9.0258x10" 

1 

-1.05068x10"^ 

2.5161xl0" 3 

2.8694xl0" 4 

-1.3673xl0" 5 

CaI 

Co 

30 

7.1551x10" 

1 

-7.2213xl0" 2 

2.4297xl0 -3 

9.08X10" 6 

~3.48xl0 -7 


. 35 

6.1612x10" 

1 

-5 .5237xl0~ 2 

3. 803xl0" 4 

2 . 8605xl0~ 4 

1.5623xl0 5 


40 

5.2087x10" 

1 

-3.5109xl0" 2 

-1, 4939xl0 -3 

3 .5115xl0 4 

-1.4808xl0" 5 
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P = corrected pC0 9 , mm Hg 

3 . * 

V = volumetric air flow rate, scfm 


pC0 2 * process air inlet pCC^, mm Hg 

A comparison of experimental data and the simulation curves for TI as a function 
of process air flow rate at 20 ASF is presented in Figure 23. 

Based on the two correlations for TI, a combined polynomial describing TI as 
a function of pC0 2 and process air flow rate at each current density can be 
written as 


TI 


S, .P + 

1,1 a 


2 

S- .P + 
2,1 a 




a 


Where 

j = 1 at 10 ASF, 2 at 15 ASF, ... and 7 at 40 ASF 

Cell Voltage 

As predicted by theory, experimental results show cell voltage was to be a 
function of current density and temperature. Data obtained for cell voltage 
as a function of current density (at 78F) were fitted to the general voltage 
current relationship as presented in the theoretical cell voltage discussion. 

Cell voltage data obtained as a function of temperature were used to correct 
the voltage current correlation on temperature values other than 78F.^ The 
temperature correction was found to be a constant 5 mv/F. The resulting 
empirical voltage correlation is 

E - 0.729 - 0.22 lnl + 0.008 I + 0.005 (T-78) 


Where 


E = cell voltage, volt 
I = cell current, amp 
T = cell temperature, F 

Figures 24 and 25 show the simulation curves and experimental data for cell 
voltage as a function of current density and temperature, respectively. 

Moisture Tolerance 

For proper EDC performance the cell must remain in moisture balance as pre- 
viously defined. The computer simulation model must, therefore, insure that 
the operating conditions are within the cell moisture tolerance. In the 
actual operation of an EDC, the cell (air out) temperature is controlled at 
a set differential temperature above the process air inlet dew point temperature. 
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FIGURE 23 TI AS A FUNCTION OF AIR FLOW RATE AND PC0 2 
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FIGURE 24 CELL VOLTAGE AS A FUNCTION OF CURRENT DENSITY 
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FIGURE 25 CELL VOLTAGE AS A FUNCTION OF TEMPERATURE 
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The temperature difference is controlled at a value determined by the nominal 
range in process air humidity conditions and current density (H 2 0 production) . 

Once the magnitude of the temperature differential is set, the cell temperature 
control scheme will maintain moisture balance based on the incoming process 
air humidity conditions. Should humidity conditions fall outside the nominal 
range, the cell temperature control scheme may be unable to maintain cell 
moisture in balance. 

The moisture tolerance of a cell initially charged with 61.5 ± 0.5% w/w Cs^CO^ 
at current densities between 10 and 40 ASF has been determined experimentally 
to be given by the following set of conditions. 

1. 9.5 ST . DPj -13 and 9.5 -1 2 - DP 2 ^13 

2. 4-Tj - DPj *14, 7 ST 2 - DP 2 *19, 9.5 - (1^ - DP X + Tj - DP 2 )/2 
and either 

13 *T 2 - DP 2 
or 

10 - (T 1 - DPj + T 2 - DP 2 )/2 

Where 

T » process air inlet temperature, F 
DP^ - process air inlet dew point, F 
T 2 = process air outlet temperature, F 
DP 2 - process air outlet dew point, F 

One of the two conditions must be met to maintain moisture balance. These 
conditions are shown graphically in Figure 26. 


Heat Removal 


For the computer simulation, it is assumed that there is sufficient heat removal 
capacity in the cooling air (flow rate and temperature) to maintain cell tem- 
perature (and, hence, moisture balance) at the required level. For a given 
cooling air temperature, the cooling air flow rate required to maintain cell 
temperature must be calculated. Two correlations are required for this calcu- 
lation: the convective heat transfer coefficient and the cooling fin efficiency. 

The correlations used were obtained from data presented in literature. ^ J 
For the geometry of the cell cooling fins attached to the anode current collector, 
the convective heat transfer coefficient is assumed to be a function of the cool- 
ing air flow rate only. The following correlation was used for the convective 
heat transfer coefficient in a rectangular channel. 
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FIGURE 26 EDC MOISTURE TOLERANCE 
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If V *1, h = 1.97V 0,3333 

If 1<V<2 .9883, h = 1.21 + 0,34 + 0.42V 2 

If 2.9883 -V -4, h = 2V 

If 4<V, h = 2.639016V 0,8 

Where 

V - cooling air flow rate per cell, scfm 

2 

h = convective heat transfer coefficient, Btu/hr-F-ft 
The following formula was used to calculate the fin efficiency 
If V<0. 32194, n= 1 

If 0.32194 <V<24 . 72 , n= 1.0744 - 0.1135V°’ 5 
If 24.72 <V, n= 2.486V -0,5 

COMPUTER SIMULATION MODEL 

The computer simulation model uses the empirical correlations obtained for 
TI, cell voltage, moisture tolerance, and heat removal to describe EDC perform- 
ance and solve EDC heat and mass balance equations as a function of 14 input 
parameters. Figure 27 is a block diagram of an EDC showing the six process 
streams entering and exiting an EDC module. The 14 input variables, or para.- 
meters, describe the inlet process streams and module-operating parameters. 

The ranges in the input variables over which the simulation model can predict 
accurately EDC performance are given in Table 5. 

Model Analysis 

The simulation model completes the description of the process air and FL 
streams by calculating component flow rates and partial pressures. 'The vapor 
pressure and relative humidity of each stream is calculated from the dew point 
temperature using a correlation obtained from the literature. ^ Cell tem- 
perature is calculated from the process air inlet dew point temperature and 
DELT1. The TI and cell voltage correlation were solved and the results were 
used to calculate module performance parameters. 

Cooling air flow rate is determined by calculating the temperature differential 
between the cell and cooling air temperatures required to remove the waste 
heat generated. The available temperature drop is assumed to be the log mean 
average temperature differential between the cell temperature and the cooling 
air inlet and outlet temperatures. Cooling air flow rate is then calculated 
by an iterative alogorithm such that AT required and AT available are within 
±0.001 F. 
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FIGURE 27 EDC BLOCK DIAGRAM 
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TABLE 5 INPUT VARIABLE RANGES 


Process Air Inlet 


pC°2 

P 0 

Total Pressure 
Dew Point Temperature 
Dry Bulb Temperature 
Flow Rate 


0.5-10 mm Hg 
2-4 Psia 
14.7 Psia 
41-70F 
47 -8° F 

0.2-0. 8 Scfm/Cell ^ J 


Process 1^ Inlet 
Flow Rate 

Dew Point Temperature 
Dry Bulb Temperature 
Pressure 


-1.3 x Stoichiometric 

10-40F 

65-75F 

0-5 Psig 


Module 

Number of Cells 

Current 

DELTlO) 


1-90 

2.44-9.76 Amp 
15-25F 


Cooling Air Inlet 
Temperature 


-Process Air Inlet 


(a) 7 OF and 14.7 Psia 

(b) DELT1 = (Module Temperature - Process Air Inlet Dew Point) 
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Mass balances for CCL, 0 2 , and H^O are used to calculate the process outlet 
stream properties. The H 2 dew point out is assumed to be the average of the 
process air inlet and outlet dew points. A complete list of input and output 
variables is given in Table 6. 

Computer Program 

The computer program was written in Fortran IV for a Univac 1108 Digital 
Computer. Figure 28 is the flow chart for the computer program. The program 
line numbers are referenced above each step in the flow chart. The subroutine 
flow charts for calculating TI, cell voltage, dew points, FLO vapor pressure 
and the solution of root functions are given in Figure 29, A listing of the 
computer program is given in Appendix C. 

The program provides checks to insure that 

1. there is a -1.3 x stoichiometric-in-current flow, 

2. process air humidity conditions are within cell moisture tolerance, 
and 

3. cooling air temperature is low enough and sufficient heat is 
generated to control module temperature. 

Model Results 

Emphasis was placed on the empirical correlations used in the computer program 
which resulted in the final simulation curves (shown as the solid lines) pre- 
sented in Figures 22 through 26. Several sample outputs showing all calculated 
variables and the program variable checks are presented in Appendix D. 

CONCLUSIONS 

The following conclusions were reached as a direct result of mathematical 
model activities. 

1. An EDC computer simulation model based on experimental results 
has been developed to predict EDC performance as a function of 
the input operating parameters for LSI's baseline cell 
configuration, 

2. The simulation model can be used for integration studies when inter- 
facing the EDC with other regenerative life support subsystems. 

3. The computer model can be used to optimize performance as a 
function of the input conditions for the cell design modeled. 

4. TI can be assumed to be a function of process air inlet pCC^, 
current density and process air flow rate over normal EDC oper- 
ating ranges. 
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TABLE 6 MATH MODEL INPUT AND 
OUTPUT VARIABLES 



Process Air Inlet: 


Inlet pCO^, mm Hg 
Inlet p0 2 , Psia 
Pressure, Psia 

Inlet Dry Bulb Temperature, F 
Inlet Dew Point Temperature, F 
Inlet Flow Rate, Scfm 

Inlet Flow Rate, Slpm 
Inlet Dry Bulb Temperature, F 
Inlet Dew Point, F 
Pressure, Psig 

Number of Cells 
Current, Amps 
DELT1 , F 

Inlet Temperature, F 


Temperature , F 
Cell Voltage, Volt 
Stack Voltage, Volt 
Current Density, Asf 
Power Output, Watt 
Heat Load, Btu/Hr 
0 ^ Consumption, Lb/Hr 
Consumption, Lb/Hr 
H-0 Production, Lb/Hr 
Transfer Index, Lb C0 2 /Lb 0_ 
Transfer Efficiency, % * 

CO Removal Rate, Lb/Hr 
Preferred Value of DELT1, F 

pC0 2 , mm Hg 
Total Pressure, mm Hg 
pH^O, mm Hg 
pN_, mm Hg 

HO Flow Rate, Lb/Hr 
C0 2 Flow Rate, Lb/Hr 
Relative Humidity, % 

0 Flow Rate, Scfm 
H 2 0 Flow Rate, Scfm 
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Table 6 


- Continued 

Process Air Outlet: Total Pressure, Iran Hg 

pC0_ , mm Hg 
pH 2 6, mm Hg 
Total Flow Rate, Scfm 
C0 2 Flow Rate, Scfm 
HO Flow Rate, Scfm 
N 2 Flow Rate, Scfm 
of Flow Rate, Scfm 
C0 2 Flow Rate, Lb/Hr 
H 2 0 Flow Rate, Lb/Hr 
Dew Point Temperature, F 
Relative Humidity, % 

Process Inlet: Total Flow Rate, Scfm 

2 pH 2 0, mm Hg 

Total Pressure, mm Hg 
pH 2 , mm Hg 
H 0 Flow Rate, Scfm 
H 2 Flow Rate, Scfm 
H 2 Flow Rate, Lb/Hr 


Anode Gas Outlet: Dew Point Temperature, F 

pH-0, mm Hg 

CO 2 Flow Rate, Lb/Hr 

H 2 Flow Rate, Scfm 

C0 2 Flow Rate, Scfm 

Total Pressure, mm Hg 

H-0 Flow Rate, Scfm 

H 2 0 Flow Rate, Lb/Hr 

Total Flow Rate, Scfm 

Total Flow Rate, Slpm 

H 2 0 Flow Rate, Slpm 

CO- Flow Rate, Slpm 

H- 2 Flow Rate, Lb/ Hr 

Total Flow Rate, Lb/Hr 

C0--H- Weight Ratio, Lb-CO^L^ H 2 

H 2 -C0 2 Volume Ratio, Ft*H 2 7Ft C0 2 

Process Air Heat Pick-up, Btu/Hr 
Process H 2 Heat Pick-up, Btu/Hr 
Cooling Air Heat Pick-up, Btu/Hr 

Flow Rate, Scfm 
Inlet Temperature , F 
Outlet Temperature, F 

Fin Convective Heat .Transfer Coefficient, 
Btu Hr Ft F 

Fin Efficiency, Fraction of Unity 


Heat Balance: 


Cooling Air: 
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5. Series flow of H 2 to a 15-cell stack does not influence transfer 
efficiency or cell voltage for anode gas pCO^ levels <640 mm Hg. 

6. The formation factor for the cell electrolyte-filled matrix was 
found to be 6.0 ±0.2. 

7. The mathematical model computer program should be expanded to 
include all components of an Electrochemical C0 2 Collection Sub- 
system and subsystem control schemes to enable total subsystem 
characterization and integration (by simulation model) with other 
spacecraft life support subsystems. 
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APPENDIX A FORMATION FACTOR CALCULATIONS 

The formation factor is defined by the equation 

f = R JR 
v' o 

Where 

f = formation factor 

R y « resistivity of electrolyte in the matrix, ohm-cm 
R q = resistivity of bulk electrolyte, ohm-cm 
The equation relating the voltage drop caused by internal cell resistance is 

e ir= ( R v )U)a) 

Where 


= internal resistance voltage, volt 

1 = compressed matrix thickness, cm 

2 

i = cell current density. A/cm 
For 55% CS 2 CO 2 charge concentration 
Ejr = 0.037 volt 

i = 0.0215 A/ cm 2 (20 ASF) 


R =4.65 ohm-cm (at 66.5 F) 

0 

1 = 0.0635 cm (0.025 in) 
Therefore 


R v = E IR /(l)(i) 

R y = 0.037/(0.0635) (0.0215) 
R v = 27.1 ohm-cm 
f = 27.1/4.65 
f = 5.8 


For 61.5% CS 2 CO 2 charge concentration 
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E jr = 0.040 volt 
i = 0.0215 A/ cm 2 

R =4.69 ohm-cm 

0 

1 = 0.0635 cm 
Therefore 

R v = 0.040/ (0.0635) (0.0215) 
Ry = 29.3 
f = 29.3/4.69 
f = 6.2 
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APPENDIX B AIR CAVITY MASS TRANSFER COEFFICIENT DERIVATION 

The air cavity convective mass transfer coefficient, k , for the EDC can be 
estimated as a function of the Reynold's number and the Schmidt number. The 
equation for this functional relationship is 



a, b and c = constants 

D. n = diffusion coefficient of C0„ in air 
AB 2 

L = diffusion length 
Re = Reynold's number 
Sc = Schmidt number 


The following assumptions can be made: 

1. Over the typical EDC operating temperature range (60-80F), the 
physical properties of air and CO^ and the diffusion coefficient 
are constant. 


2. The Schmidt number is approximately 1.0 (<5% error). 

3. The Reynolds number is a function of the bulk air flow rate 
only. 

Using these assumptions the equation for k reduces to 

c 

k c = A v b (1) 

Where 


A, b - constants 

V = volumetric air flow rate 

The value for k can be found experimentally by measuring the C0 2 removal 
efficiency audience, removal rate as a function of air flow rate for a low 
inlet pC0 2 level (-2.0 mm Hg where air cavity transfer is rate-limiting). 

As previously shown, the equation for CO^ removal is 
n co 2 - Ck c ) CA) CM C o 2 ) C P co 2 - pco 2 ° )/RT 

Solving for k yields 
c 
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\ - N C0, RT/AM C0, ( P C0 2 - P C0 2°> • C2> 

In the air cavity mass transfer limiting region, pCCL 0 is assumed to be zero. 

The value of k £ can then be calculated by experimentally measuring N fQ for 
a given inlet pC0 2 and air flow rate. The values for k at each air 2 flow rate 
tested (for a given inlet pC0 2 ) were least squares fit £o equation (1) to 
give the value for A and b. 
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APPENDIX C 


COMPUTER PROGRAM DOCUMENTATION 


C-l 




Chi/US 16 OCT 73 


WU; 2 


I. IFt SYSTEMS U34900.0) 


21 S 40! 2? 


DATE % 16 OCT 73 


AH RUE 

FOkIkAN. V - SEPTEMbER 20 r 1973 TUESDAY » OCTOBER 16 . 1973 AT 21 140: 24 

HIbhLSr S«E • AJDHESS: 03110 


UOulul 1* KEALlV.I 

00ulU3 2* COMMON A . 8 . C . r r Ft r F » G . CPA . OHC » N r T2 » T5 » T6 . HC . ETA » S l 5 . 7 ) 

0 0 0 i U 4 3. DATA S/ 1 . 75187 * - . 51 1449 . . 0 7760 73 . “ . 005925*4 . .000179862* 

0001U4 4. SI *42306. -.271^35. . 0210 187 . 0 0 0 34839 . -1 . 8475E-5 * 

OUo 1 04 5. #1. 00190 1 1062 1 00814(i8 » * 0 0 196359 * -a. 5832E-5 » 

000104 t>. S.9o2Sa»“.105uftt» .0025161 » . 00 028694 . -1 . 3673E-5 » 

00U1U4 7. ft .71551 »-. 07221 3 * » 0024297 r 9. 08E-6. -3.48E-7# 

0 0 0 1 U 4 a. 61 6 12 .-. 055237 . . 0003803. . U 0 028605 » -1 . 5623E-5 . 

000104 9. «. 52087 03bl 09 »-*• 00 14939. . 0 00361 1 5 » -1 .48P8E-S/ 

000 1U4 lu. 4. A.B»C»D»EE»F»b»CPA.CPH/3.2437814» . 00566826 * 1 . 170P379E-8 . 

UOU1U4 11. S. 0021878462. 5. 2 19603 *273. 16. 2. 3025851 » 1.0765. 1.0666/ 

uOOl 17 12. I IN = 5 

00U120 13. IOUT = 6 

00O120 14, C WHEN NELAG=1 AND CERTAIN PARAMETERS ARE OUT OF RANGE. 

U0ijl20 15. C ONLY THE ASSOCIATED MESSAGES AND THE INPUT ARF PRINTED- 

U 00 120 16. C WHEN NFLAG=0 A COMPLETE PRINTING IS DONE. 

I 000121 17. 1 REaD(IIN.7O0i.enD = 9999) PC 1 . PO IPSA . P1PSA » T1 # DPI * V 1 . 

M 000121 18. SV3SL.T3.DP3.P3PS6.N. I. DFLT1 » T5 . MFLAG 

u 0 l.j 1 42 19. 7001 FORMAT ( 6F10 .2/4F10 • 2/4F1 U. 2 » 1 10 > 

000143 20. WRITE { IOUT. 230 

000145 21. 230 FORMAT ( »!**♦*****♦* INPUT DATA********* * * ) 

000146 22. WRITE! IOUT » lOn > PCI. P01PSA. PI PSA. Tl. DPI .VI 

000156 23, 100 FORMAT! »0CARD 1* PROCESS AIR IN! 1 /* PCI P01PSA*. 

00U156 24. «* P1PSA T 1 DPI V1*/6FJ0.2> 

OUu 157 25. WRITEdOUTMO?) V3SL» T3 » DP3 .P3PSG 

000165 26. 102 FORMAT! ♦OCARU 2 : H2 IN:V» V3SL T3 PP3'# 

000165 27. 4* P3PSG*/ 4F10.3) 

UOu lob 28. WRITE ( IOUT. 104) N. I » CELT 1 . T5*NFtnG 

OOoi 75 29. 104 FORMAT (*0C AND 3 ! MODULE S * / * N I DFLTl 1 » 

000175 30. *• TS NFLAG*/F10c0,3F10.2.I1u) 

0 0 0 1 7 6 31. WRITE ( IOUT. 220 ) 

0 0 0 2 u 0 32. 220 FORMAT! » u***********RESULTS************ * ) 

U002U1 33. NErKOR=0 

00u2Ol 34, C 

uOOeiol 35. C INLET PRoCf SS wIK CONVERSIONS: 

Oli02ul 3o. C 

001-202 37. Pl = PiPSA*51.70ii7 

000203 38. PQi=:POlPSA*bl . 7007 

U002u 4 39. PWi=PHTO!DPl) 

OOGcUS 40. PNl=Fl~POl“Pvvl-PCl 

000*06 41. VW1=PW1/P1*V1 

continued- 
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Figure C-l continued- 


00020? 

42. 


VtU-POi/Pl + vl 

UU0210 

43. 


F\n1 = VW1*10.U1*oO./386.7 

UOu^il 

44. 


FC1=V1 *PC1/H 1*44. 01 *60. /38b. 7 

00O212 

4b. 


VC1=V1*PC1/P1 

OU0213 

4b . 


RH i— PW i /PHTO ( T 1 > *10 0 . 

uaL .^13 

47. 

C 


U002i3 

4b. 

C 

INU5T H2 CONVFkSIONS: 

ouu2i3 

49. 

L 


000214 

50. 


V3-V5SL/20.32 

0Uo<r.l5 

51. 


PifJ5=PHT0(DP3) 

UOO^ib 

52. 


P3=(P3P56+14.7> *5l.7UU 7 

000217 

55. 


PH3=P3-PW3 

000220 

54. 


VW3=PW3/P3+V3 

000221 

55. 


Vh5=V3-Vw3 

000222 

56 . 


FH3=VH3*60 • *2 • 01 b/306.7 

(W 02.22 

57. 

c 


000222 

Sti. 

c 

MODULE PaR ^ METER CALCULATIONS: 

UOl i2<;2 

59. 

c 


00 u22 3 

60 * 


T2=(jpi+neLTi 

U00224 

61. 


E=.729-.22*Al 0. 3 ( I)+.OO8*I+.OOb*tT2-70. ) 

000225 

62 . 


ESIACK=N*E 

oon^fe 

65 • 


CUNUEN-I/ .?44 

000221 

64 . 


POwEK=£STACk*I 

000250 

65. 


HEATLO=N*I*( 1.25-E)*3.419 

U0u23l 

66 • 


OCON=6.5603t>4*M*I 

000232 

67. 


HC0N=0C0im*2.016/16. 

UU0233 

6ti, 


WPnOu = OCON + PICON 

000234 

09. 


TI-T1C0R (PC1»V1/N*CUK0EN) 

0002 Ob 

70. 


TE=TI/ . 0275 

000206 

71. 


CTKanS=TI*OCOM 

0 0 0 2 0 7 

72. 


IF l Fh3 . PE • 1 » 3* t-iCON ) GO TO 1213 

U00241 

73. 


V 3PSL= I. 3*HCQN+V5*386. 7*28.32/ lvH3*2 #016*20.) 

000242 

74. 


WP1TE( IOUT»213)V3MSL 

000245 

75. 

213 

FORMAT ( • 0 INSUFF IC IENT HYDROGEN. '»F10.2r* SLPM IS PEPUIRFO.M 

UUl/246 

76. 


NEHR0R=1 

00U246 

77. 

t 


000246 

7tt • 

C 

COOLING AIR AA'U HEAT RALANCF CALCULATIONS: 

000246 

79. 

c 

INC* ETA AND To ARE CALCULATED IN FUNCTION ERT.) 

00o24b 

00. 

c 


000247 

01. 

1213 

DHA=V1*CPA*<T?-T1) 

000260 

82. 


DHH=V3*CPH* (T2-T3) 

U0l*25l 

85. 


OHC=HF«TLD-DHA-DhH 

000252 

04. 


I F IDHC • GE*0 • ) GO TO 098 

000254 

65 . 


WRITE l I OUT* 789>PHC 

U0u2b7 

06* 

789 

FORMAT (•ODHC " ' » IPG 10 • 3 * * HTU/HR.V INSUFFICIENT HEAT •» 

000257 

87. 


fit IS GENERA TFO 10 MAINTAIN MODULE TEMPERATURE • * ) 

U0u2oU 

0d. 

676 

V5=u. 

U 0 0 2 o 1 

89. 


HC — 0 . 

U0u2o2 

90. 


T6-0. 

U0O2b3 

91. 


ETa=0. 

0002O4 

92. 


nerror-i 

0 0 0 2 6 6 

95. 


GO To 1430 

00L2ob 

94. 

696 

CONTINUE 

U 0 0 2 6 7 

95. 


lF(T2-T5.GE..37b*DHC/N>G0 TO 765 
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o0u2'U 
000273 
000274 
u0ii«i75 
0002/6 
U 0o277 
000300 
UUO 3ol 
OtH.oul 
00(j 3u2 
U0u3U3 
0003U3 
O0U3U3 
00GOU3 
U0U3U6 
UUUJ07 
U0031U 
001,311 
U0o312 
00u313 
U00314 
U00315 
000316 
000317 
0Ui320 
000321 
U0U322 
000323 
000324 
00 1.324 
000324 
000324 
000325 
00 0326 
U00327 
000330 
000331 
U0U332 
000333 
U0U334 
00U335 
000336 
000337 
000340 
000341 
000342 
000*344 
000345 
000345 
000345 
000345 
00U346 
00034b 
UUo3b0 


9o. 


ttP I Tfc. ( I OUT * 967 ) 

97. 

967 

FORI*AT( *0 INSUFFICIENTLY LOW cooling AIK TEMPERATURE 

9b. 


GO TO 67b 

99. 

76b 

CONTINUE 

100. 


VA = 1.01*UHC/<CPA*N'MT2“T5) ) 

101. 


VH=2. *VA 

102. 


V5=n*ROOT(VA»VH»EPT»0. » .ij 05» .flQlrlU) 

103. 

1430 

CONTINUE 

104. 

C 

trial VALUE OF HP4 ; 

105. 


UP4-T2-1 1 • 

10b. 


DO 937 NNN*J=1»2 

10 7. 

c 


106. 

c 

OUTLET H2 calculations: 

109. 

c 


110. 


PW4=PHT0(DP4) 

111. 


FH4=FH3-hC0N 

112. 


FC4=CTkANS 

113. 


VC4=FC4*386. 7/44 .0 1/60 . 

114. 


VH4=VH3-HC0N*366.7/6G./2.016 

115. 


P4=P3 

116. 


V W 4— P W 4 / ( P 4 “ p i« 4 ) * ( VC4 + VH4) 

117. 


FW4 = \/W4*bO* *18.01/386. 7 

116. 


V 4=VH4+VW4+VC4 

119. 


CHWRT0=FC4/FH4 

120. 


HCVRTO=VH4/VC4 

121. 


VH4SL=VH4*28,32 

122. 


V4SL=V4*28.32 

123. 


VC45L-VC 4*26 * 32 

124. 


F4=FC4+Fw4+FH4 

125. 

c 


12b. 

c 

OUTLET PROCESS AIR CALCULATIONS: 

12 7. 

c 


126. 


P2=P1 

129. 


FC2=FCl-CT«ANS 

130. 


VC2=FC2*386.7/b0./44.01 

131. 


V02=V01-OCGi j*366 *7/32 . /60 . 

132. 


VN2=PN1/P1*V1 

133. 


VW2=VWl+VW3+WFKOU*38b.7/6O./l0.Ul - VW4 

134. 


V2=V02+VN2+VW?+VC2 

135. 


PW2=P?*VW27V2 

136. 


PC2=P2*VC27V2 

137. 


DP2=DEwT(PW?) 

136. 


RH2=PW2/PHTQ(T2)*100. 

139. 


FW2=VW2*ia* 01*60-/386.7 

140. 


OP4=(npi+oP2J#,5 

141. 

937 

CONTINUE * 

142. 


0ELT1P=?2.-T1+JP2 

143. 


T2P=DP1+DEUT1P 

144. 

C 


145. 

C 

moisture balance checks: 

146. 

C 


147. 


IF(ABS(T1-0P1-11 .25) .LE. 1. 75 

14B. 


A. AND. AHSIT2-DP2-11. 25) .LE.l .75)00 TO 30 

149. 


IFITL-DP1.LT.a,)GO TO 1210 
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0 

1 

tn 


000432 

150. 


UU0334 

131. 


UOOObb 

152. 


00u3o0 

133. 


U0G3ol 

154. 


00ii3ol 

15b. 


uQD3o3 

156. 


00i)3b5 

157. 

654 

OOtiOob 

158. 

28 

UOu 3?0 

159. 

543 

UOOJ/I 

16U. 

29 

UQ0372 

161. 

30 

000373 

Ib2. 


000373 

163. 

C 

u0u373 

164. 

C 

00(j373 

165. 

c 

U0O376 

166. 


00u37b 

167. 


00U413 

166. 

106 

000413 

169. 


00U413 

170, 


U0U414 

171, 


OOU427 

172 , 

101 

000427 

173. 


0QU427 

174. 


000430 

175. 


000430 

176. 


000446 

177, 

111 

000446 

178. 


000446 

179. 


000447 

180. 


000460 

181. 

103 

UGO460 

182. 


00U4ol 

183. 


0 0 U 4t>l 

164 « 


U00503 

183. 

110 

O0G3u3 

186. 


U0(J bu3 

187. 


O0ubu3 

186. 


OOU304 

189. 


UUUbll 

19U, 

109 

OOubll 

191. 


UOOD12 

192, 


000321 

193. 

107 

000521 

194. 


000322 

195. 


000330 

196. 

4 56 

UOlJ 330 

197. 


000330 

198. 


OOubbl 

199. 


000332 

200. 

1210 

000334 

201. 

432 

U0U335 

202. 


00UbO6 

203. 

1211 


lF(Tl-uPl.ftr.m.)GO TO 1210 
IF!T2-OP2,LT,7. )G0 TO 1211 
IFIT2-DP2.GT.19- IGO TO 1211 
DPTCJAV = .5*!T1-DP1+T2-0P2> 

I F ! AHS ( 0PT0A V— 1 1 • ) • LE • 1 • b 

& .AND • (T2-DP2 .G£. 13. .OR. 0PTPAV .GE. 1 -i . > J RO TO 30 
WRITE ( IOUT » b5A ) 

FORMAT ( * 0 AVERAGE PROCESS AIR OEw POINT DEPRESSION OUT OF RANGE. •> 
WRITE! 10UTr543> 

FORMAT ( * OELECTROLYTE MOISTURE BALANCE IS NOT MAI NT A T NED . * ) 

nekror=i 

CONTINUE 

if(nerror*nflag.ne.o>go TO 1 
MAJOR OUTPUT satements 

WRITE ( lOuT * 1 06 ) T2 » E * ESTACK * CURQFN * POWER * HEATLO » OCON * HCON » 

«WPROl)» Tl* TE* CTRANS 

FORMAT ( «OMOOULg: */BX» *T2* *9X* « E ' * 4X * 1 ESTACK * *4X» *CURDEN» * 

&5X » * POWER 1 *4X» ♦H£ATL0»/6F1Q,3//6X» ‘OCON* *6X* *HCON» »5X* »WPROD* » 
S0X* »TI ■ * 8X* * TE’ *4X» »CTRANS , /3F10.5»3F10.4) 

WRITE! IOUT * 101>P01*Pl»PWl»PMi»Fwl*FCl»HHl*V01»VWl 

formati’oprocfss air in: */* poi pi pwi*» 

&7X* »PN1 ’ *7X. *Fwl* »7X* *FC1*/6F10,3//7X**RH1 * ♦ 7X » * VOl * r 7X r * VW1 . * / 
83F10.3) 

WRITE C IOUT » 111 )P2*PC2»PW2» V2* VC?* VW2» VN?» V02*FC2»FW2» 

SDP2*RH2 

FORMAT t * OPROCFSS AIR OUT l ♦ / » P2 PC2 PW2 ' * 

S0X» 1 */2* »7X* *VC2* »7X**VW2*/4F10.2*2F10.4//7X* *VN2* *7Xr *V02** 

&7X » ' FC2 * * 7X * » Fv«2 * * 7X * *QP2* *7X» » RH2 ' /6F10 . 3 ) 

WRITE! IOUT* 103 )V3*Pw3* P3*PH3 »vw 3* VH3*Fh3 

FORMAT! *0H2 INI •/* V3 PW3 P3 PH3* * 

«» VW3 VH3 FH3»/F10.4*3F10.2»F10.6*F10.4»F10.6) 

WRITE! IOUT* llO)nP4*PW4»FC4* VH4* VC4rP4» VW4*FW4, V4* V4SL* 

& VH4SL t VC4SL * FH4 * F4 # CHWRTo » HCVRTO 
FORMAT! *OANODE GAS OUT:*/* DP4 PW4 FC4'» 

8 7X* *VH4* *7X* »VC4* *8X* * P4 * /3F 10 . 3 • 2F1 0 • 6 » FI 0 . 2//7X * »VW4» *7X*»FW4» » 
8ttXr *V4* »6X* 'V4SL' »5X» *VH4SL* *5X» * VC4SL * /2F1 0 . 6 * 4F1 0 . 4// 

«7X» * FH4 * »8Xr »F4» »4Xr * CHWRTO* »4X* * HCVRTO ’ /F 10 . 6 * 3F1 0 . 4 ) 

WRITE! IOUT * lUQiDHArOHHrDHC 


FORMAT! *UHEAT BALANCE:*/* 
«3Fin.3) 

DHA 

DHH 

DHC » / 

WRITE! IOUT *iU7)V5*T5*T6*HC» ETA 
FORMAT ( * UCOOLING AIr:»/« 

V5 

T5» # 



R * Tb HC ETW5F10.3) 

WRITE! IOUT *456) PELT IP »T2Pf DPI » BELT IP 
FORMAT ! * URREFEHPED VAUJF OF DELTl IS»*F6.P* 

H* UNDER THESE CONDITIONS*/* SO THAT » »FR. 2 » 

8* = T2 = DPI + PELT 1 =»*F6.2*» +»»F6. ?**.») 

GOTOI 

WRIT£(I0UT*432) *■ 

FORMAT! *0INLET PR0CES5 AIR DEW POINT DEPRESSION OUT OF RANGE.*) 

GO TO 28 

WRITE t IOUT » 321 ) 



Figure C-l continued- 


000540 

204. 

321 

FOkMAT( 'OOUTLFT process A IR DEW point depression out OF range. *> 

U00541 

20b. 


GO TO 2 H 

000342 

20b * 

99S9 

STOP 

00o342 

207. 

t 


00{;542 

2U6. 

C 

WATER VAPOR PRESSURE (PHTO) IS A FUNCTION OF PEW POINT (TF> 

006342 

200. 

C 


U0ijb43 

210. 


FUNCT IONPHTO t TF ) 

OG0546 

211. 


REAL N r I 

00lib47 

212* 


COMMON A »B .C »P» FE.F.G. CPA. DHC » N» T2» T5» T6» HC » ETA 

UOUbbO 

213. 


TC=b.*(TF-32, )/9. 

UOubbl 

2X4. 


X— 374 . X 1-TC 

UGlibb2 

21b. 


PHT0=EXP(G*(EF-<X/(TC+F) )*( A+R*X+C*Y**3>/< 1.0+D*X) ) ) 

U0libo3 

216. 


RETURN 

UQ0bb3 

217. 

C 


000353 

2l«. 

C 

UEW POINT ( OEWT ) IS A FUNCTION OF WAT^R VAPOR PRESSURE CP)* 

00ubb3 

219 . 

c 


UUobb4 

220 . 


FUNCTION DEwT(P) 

U00bb7 

22 1 . 


REAL n. I 

OOoboO 

222 . 


COMMON A »6»C . P . Ft » F # G. CPA »UHC »N» T2 » Tb » 76 . HC ♦ ETA 

OOObbl 

223. 


EXTERNAL ROOT. PHTO 

000562 

224. 


XX=ALOG(P) 

000563 

22b. 


TA=-2.4+20.25*XX+1.522*XX^*2 

000564 

22b. 


TB=TA+.l 

00ubfe5 

227. 


OElwTZROOT ( TA » TH » PHTO » P . .005# .005*5) 

UOObbb • 

226. 


RETURN 

UOObbb 

229. 

c 


oOFbbb 

230. 

c 

TRANSFER INDEX (TICOH) IS A FUNCTION OF INLET PC02 (P)» 

000566 

231. 

c 

AIRFLOW PER CFlL <AF)» AND CURRENT DENSITY <Cn). 

000567 

232. 


FUNCTION TICOP(P.AF.CD) 

000572 

233. 


REAL N. I 

00ob73 

234. 


COMMON A . R# C .P»EE* F » G» CPA* DHC .N. T2» Tb. T6 » HC » ETA . S ( 5 » 7 ) 

U0U574 

23b. 


p A— P* ( AF/ *44 ) **( • 45** • 0 (J4b*p*P ) 

00ub75 

23b. 


J=C0/5.-l. 

OOGb 76 

237. 


IF(J.LT.I) J=1 

UOQoOO 

23b. 


TI1=PA*(S( 1. J)+PA*(S<2. J>+PA*(S(3» J)+PA*tq(4» j)+PA*S(5.ji ) > > > 

OOCiOjl 

239. 


IF ( J • LT « 7 ) GO TO b 

0Q06U3 

240. 


TIC0K=TI1 

00U6U4 

241. 


RETURN 

OOObOS 

242. 

a 

J=J+1 

GOGoUb 

243, 


TI2=PA*<S(1» J)+PA*tSC2. J)+PA*(S13» J)+PA*(^(4.J)4PA*S(5.J) ) > ) ) 

0006U7 

244. 


A 1 1 — J 

UOU610 

24b. 


TEHP-CD/5 .“A 1 1 

QOUol 1 

24b. 


TIC0R=TI2*TERP+TI1* l 1 .-T£RP) 

000612 

247, 


RETURN 

00O612 

246. 

c 


000612 

249. 

c 

FUNCTION ROOT RETURNS THE NUMBER “ROOT- SUCH THAT 

0006X2 

2b0. 

c 

YFCN l ROOT ) = W + 0R- YTOL » AND SUCH THAT THE LAST 

000O12 

251, 

c 

CHANGE IN -ROOT- IS LESS THAN XTOL. MAXIMUM NUMBER 

000612 

252. 

c 

OF ITERATIONS - K. A SECANT FOLLOWING METHOD IS 

U0uoi2 

253. 

c 

USED STARTING wITH X-VALUES XI AND X2. 

0Q(jol2 

254. 

c 


000013 

255. 


FUNCTION ROOT (Xl » X2 » YFCN . W . XTOL . YTOL . K ) 

00 06X6 

256. 


EXTERNAL YFCN 

U0UOX7 

257, 


X A- X 1 



Figure C-l continued- 


00(ib20 

256. 


XP-X2 

001.621 

259. 


FA=YFCm(XA)-W 

O00o22 

260. 


FR-YFCn ( XU ) - w 

UU0o23 

261. 


DO 91 I POOT-1 r k 

UOt'6^6 

262. 


XN— (F A*XB**FB*y a ) / (FA-FB ) 

UUUt>27 

263. 


FN- YFCN < XN)-W 

U0U»$U 

264. 


1F(A6S(FN) .LT.tTOL .and. 

00o650 

26b. 


«ARS<XN-XH) .LT.XTOUGO TO 99 

i)00bo2 

266 • 


IF < IKOOT.GE.K) «6lTEtI0UT#200)XlfX2rW»X TOL t Y TOL » IROQT 1 

00u632. 

267. 


H X A f X 6 » >1 N t F A r F P » F N 

OOObbl 

266. 

200 

FORMAT t ♦ONONCOfiVERGENCE IN ROOT' » 1P3E1 2 . 4/2E 1 ? . 4 r I 3/bF 12 . 4 ) 

00(.ibb2 

269. 


XA=Xri 

00u6b3 

270. 


XB = XN 

U0D6b4 

271. 


FA = F6 

OUObbS 

272. 


FR=Fn 

OOlJObfc 

273. 

91 

COnTINUF 

OOObbO 

274. 

99 

r00T=Xn 

UOlbbl 

27b. 


return 

OOObbl 

27b. 

C 


UOi-Gbl 

277. 

C 

function ert calculates cooling air exit temper aturf (T6> 

000661 

276. 

C 

FROM THE ENTRANCE TEMPERATURE (T5>» HEAT DUTY < DHC ) t 

OOubol 

279. 

C 

ANN COOLING A IK FLOW PER CELL (tf). FROM T2* T5» AND T6 

U0li6bl 

260. 

C 

THE LOG MEAN MODULE-COOLING AIR TEMPERATURE DPlVING 

UOl/bbl 

261. 

t 

FORCE IS CALCULATED. THE DRIVING FORCE REQUIRED TO 

UObbol 

262. 

C 

TRANSPORT the hfat through the fins AND INTO THE 

U OLI t>61 

263 • 

C 

COOLING AIR IS CALCULATED. THE DISCREPANCY RFTWEEN 

GOObbl 

284. 

c 

THE AVAILABLE AMD REQUIRED DRIVING FORCES IS RETURNED 

OOL'bbl 

265, 

c 

THROUGH THE VARIABLE ERT. FIN EFFICIENCY (ETA) AND 

UOUODl 

286. 

c 

CONVECTIVE HEAT TRANSFER COEFFICIENT (HC) ARE ALSO CALCULATED. 

00LO61 

28/ • 

c 


0006b2 

288. 


FUNCTION ERT ( V ) 

00( ! 6b5 

289. 


REAL N»I 

UORbbb 

290. 


COMMON A»8»CrDrEE»F»GrCPttFDHC»N*T2»T5»T6fHC»ETA 

000bb7 

291. 


IF ( V.LE. i« ) HC=1 .97+V** ( .3333) 

0006/1 

292. 


JF ( V .GT * 1 • .AND, V.LT.2.9863)HC = 1,21 + V*( ,34+.42*V) 

000673 

293. 


IFIV.GE. 2*9883 .AND. V .LE « 4 • ) HC-2 • *V 

0006/5 

294. 


IF(V.GT.4. )HC=2.b390l6*V**( .8) 

000677 

29b. 


ETA=1. 

00u7u0 

296. 


IF(V.GT. .32194 .AND. V .LT . 24 . 72 > E T A=1 . 1)644- . U35*S0RT l V ) 

00u7u2 

297. 


iFtV.Gh .24.72) ETA=2. 486/SORT (V) 

000 /U4 

298. 


T6 = T5+DHC/U*CPA*N') 

0007u5 

299. 


DTKEQ=DHC*( .376+3. 876/(HC»ETA) ) /N 

000/06 

300. 


IF (T6.LT .T2 ) GO TO 86 

000/ 10 

301. 


ert=-dtreq 

oou /11 

302. 


RETURN 

OOo /12 

303. 

66 

OTAVA=(T6-T5) /AlOGt (T2-T5>/(T2-T6) ) 

UUO 713 

304. 


FRT-DT ava-dtre • 1 

000714 

30b. 


RETURN 

OOL /lb 

30b. 


EMU 

1 ulf 

1 gnostic 

• 

r 




2.5820 SECONDS. 

END OF 

rum: 

gate: 16 

OCT 73 TIME.: ai;52:57 CPU TI'-'E » 2.59 SFC. PAGES I 6 CARDS: 


o out 
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TABLE C-l NOMENCLATURE AND UNITS 


Composite Variable Names 

First Prefix P - (Partial) Pressure, mm Hg 

V - Volume Flow Rate, Scfm 70F, 760 mm Hg 
F - Mass Flow Rate, Lb/Hr 
DP - Dew Point, F 
T - Temperature, F 
RH - Relative Humidity, % 

Second Prefix C - Carbon Dioxide (CO.) 

0 - Oxygen (0O 2 

N - Nitrogen TN 2 ) 

H - Hydrogen (H 2 ) 

W - Water (HO) 2 
(None) - Total Stream 

Stream Number 1 Process Air Inlet 

2 Process Air Outlet 

3 H 2 Inlet 

4 H 2 Outlet 

5 Cooling Air Inlet 

6 Cooling Air Outlet 

Units Suffix (None) Units of First Prefix 

PSA Psia 
PSG Psig 
SL Slpm 

Additional Variable Names 

A, B,C,D,EE,F,G 
All ,J 
CD 

CHWRTO 
CPA 
CPH 
CTRANS 
CURDEN 
DELT1 
DHA 
DHC 
DHH 
DPTDAV 
DTAVA 
DTREQ 
E 

ERT 

ESTACK 


C-8 continued - 


Constants in pH^O Dew Point Equation 
Indices Determine S Constants 
Current Density 

CO /H 2 Weight Ratio, Lb CO /Lb H 

Volumetric Specific Heat, Air, Btu/Hr-Scfm-F 

Volumetric Specific Heat, H 2 , Btu/Hr-Scfm-F 

C0 2 Transferred, Lb/Hr 2 2 

Current Density, Amp/Ft 2 

Control Variable * T2-DP1, F 

Enthalpy Gain of Inlet Process Air, Btu/Hr 

Enthalpy Gain of Cooling Air, Btu/Hr 

Enthalpy Gain of Inlet H 2 , Btu/Hr 

Process Air Dew Point Depression, Average, F 

Log Mean Available Temperature Drop, F 

Required Module - Cooling Air Temperature Drop, F 

Cell Voltage, Volt 

Difference between DTREQ and DTAVA 

Stack Voltage, Volt 



Cifc Systems, Jhc. 


Table C-l - continued 


ETA 

FA,FB,FN 

HC 

HCON 

HCVRTO 

HEATLD 

I 

N 

NERROR 


NFLAG 

OCON 

P 

PA 

POWER 

S 

TA,TB 

TC,TF 

TE 

TERP 

TI 

T2P 

Til 

T12 

VA,VB 

WPROD 

X1,X2,XA,XB,XN 


Cooling Fin Efficiency, Dimensionless, 0-1 
Corresponding Dependent Values 2 

Cooling Fin Heat Transfer Coefficient, Btu/Hr-Ft -F 
Consumption, Lb/Hr 

Hj/CO Volume Ratio, Scfm H 2 /Scfm C0 2 

Net Heat Produced, Btu/Hr 

Current, Amp 

Number of Cells 

Check Violation Indicator 

0 No Violation 

1 Violation 

Program Option Integer = 0 or 1 
0 2 Consumption, Lb/Hr 
Inlet pCO- 

Inlet pC0 2 Corrected for Effect of Air Flow 
Electrical Power Produced, Watt 

Array containing Coefficients of PA in TI Correlation 
Trial Values of Dew Point, F 
Dew Point, C, F, respectively 
Current Efficiency, % 

Current Density Interpolation Variable 
Transfer Index, Lb C0 2 /Lb 0 . 

Preferred Value of T2 

TI at one of the seven current densities just below 
the desired current density 

TI at one of the seven current densities just above 
the desired current density 
Trial Values for V5 
Water Production, Lb/Hr 

Trial Values of Independent Variable in Convergence 


NOTE: Additional variable description is present in the comments- of the source 

program listing. 


C-9 
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APPENDIX D 


COMPUTER PROGRAM SAMPLE OUTPUTS 


D-l 




Cife Systems, Jhc. 


CARD 1* PROCESS AXE IN: 


PCI 

P01PSA 

P1PSA 

T1 

DPI 

VI 

3*00 

3.10 

14.70 

57.00 

46.00 

40.00 

CARD 2: H2 

IN: 





V33L 

T3 

DP3 

P3PSG 



10.500 

70.000 

46.000 

5.000 



CARD 3: MODULE: 





N 

I 

DELT1 

T5 

NFL AO 


90. 

4.88 

20.00 

57.00 

1 


a 

"RESULTS**** 

If if If If 




MODULE: 






T2 

£ 

ESl'ACK 

CURDSN 

POWER 

HEATLD 

66.000 

0.359 

32.33B 

20.000 

157.80B 

1337.405 

OCON 

HCON 

Vi PROD 

TI 

TE 

CTRANS 

0.28901 

0*03641 

0.32542 

2.1713 

78.9571 

0.6275 

PROCESS AIR IN: 





FOl 

PI 

pin 

PN1 

FW1 

FC1 

160.272 

760.000 

7.929 

588.799 

1.166 

1.078 

RHl 

VOl 

VW1 ' 




66.609 

8.435 

0.41? 




PROCESS AIR OUT: 





F2 

PC2 

PW2 

V2 

VC 2 

VW2 

760-00 

1.26 

10.15 

39.97 

0*0660 

0.5336 

YN2 

V 02 

FC2 

FW2 

DP2 

RH2 

30.989 

6*577 

0.451 

1.491 

52.614 

62.008 

H2 IN: 






V3 

Pvi 3 

P3 

PH 3 

lift 3 

V1I3 

0.37 OR 

7.929 

1019. 

1011. 

2.886E-03 

0.3679 

ANODE GAS 

OUT: 





DP 4 

PW4 

FC4 

VH4 

VC4 

P4 

49.307 

8.971 

0 *62R 

0.251461 

0.09189? 

1018.50 

Ytf 4 

PS 4 

V4 

V4SL 

VH4SL 

VC4SL 

0.003051 

0*008527 

0.3464 

9. 8103 

7.1214 

2.6025 

FH-1 

F4 

CHWRTO 

HCVfiTO 



0 *078657 

0.7147 

7.9780 

2.7363 



MEAT BALANCE: 





DHA 

DHH 

DHC 




367.540 

-1 .582 

951 .527 






D-2 





FII3 

0.1151 
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COOLING A Ili: 


v 5 

To 

T6 

HC. 

ETA 

626.416 

37.000 

58-411 

12.461 

0.765 

PREFERRED 

VALUE OF : 

DELT1 IS 17. 

61 UNDER 

THESE CONDITIONS 

60 THAT 66 

.61 = T2 • 

= DPI + DELT1 = 46.00 

+ 17.61. 

a-#*#* 

I PUT DATA********** 



CARD 1* PROCESS AIR 

IN: 



PCI 

PCI PS A 

PI PSA 

T1 

DPI VI 

o. 00 

3.10 

14.70 

51.00 

46.00 0.80 

CARD 2: ii2 

IN: 




V35L 

T3 

DP3 

P3PSG 


0.300 

70.000 

46-000 

0.000 


CARD 3: 0 ODORS: 




H 

I 

DELil 

X5 

NFL AG 

1. 

■if - ' 1 *' m mr ic’vi *!<• ■>* vS* I't * 

2 • 44 

20.00 

bl .00 

1 

M ' KA b U I/i 1 ^ ir ^ -rt ' v_ ** _J r ■* r ■ ■ 




-6.36 BTU /HR • 



INsUFFIOH;?! 

i ( HKA i’ IS 

GENERATED 

TO MAINTAIN MODULE TEMPER A TO. 




CARD 1» PROCESS AIR 

IN: 



PCI 

P01PSA 

PI PS A 

ll 

DPI VI 

2.80 

3.10 

14.70 

66.00 

57.00 0.44 

CARD 2: H 2 

IN: 




V3SI- 

To 

DP6 

P3PSG 


0 * 0 0 

70.000 

46.000 

0.000 


CARD 3: MODULE: 




N 

I 

DELT1 

To 

NFLAG 

1. 

4.80 

20.00 

68.00 

0 

■'r*!,- #-$r ■*■■'& vr i’ 

HE S U L T S -x 



average process air 

DEtti POINT DEPRESSION 

OUT OF RANGE. 

electrolyte 

MOISTURE 

BALANCE IS 

HOT MAINTAINED. 

MODULE: 





T2 

E 

ESTACK 

CUR DEN 

POWER H3A7LD 

7? .000 

0.414 

0.414 

20.000 

2.022 13.043 

OC OK 

BOOK 

ViPiiOD 

XI 

T2 C TRANS 

0.00321 

0.00040 

0.00562 

2-0677 

75.9163 0.0067 


PROCESS AIR 
POl 

1 60 - 2? 2 

IN: 

PX 

760.000 

PV/1 

11.905 

PN1 

585.025 

FW1 

0.019 

FC1 

0.011 

RH1 
67 • 886 

VOl 

0.093 

V Vi 1 

0.007 

D-3 





JCife Systems, Jnc. 


PROCESS AIR OUT: 





F2 PC2 

PW2 

V2 

VC2 

VK2 

760*00 1*11 

14®04 

0*44 

0.0006 

0.0081 

VN2 702 

FC2 

FN2 

DP2 

RH2 

0.339 0.092 

0*004 

0.023 

61.623 

59.103 

H2 IN: 





V3 PW3 

P3 

PH3 

VW3 

VH3 FH3 

1.059E-02 7.929 

760.0 

752.1 

1.105F.-04 1. 

048E-02 3.270E-03 

ANODE CAS OUT: 





DP4 P W4 

FC4 

VH4 

VC4 

P4 

59.311 12.894 

0.007 0 

.009189 

0.000982 

760.00 

V vj 4 FV« 4 

V4 

V4SL 

VR4SL 

VC4SL 

0.000176 0.000491 

0.0103 

0.2930 

0.2602 

0.0278 

FH4 ?4 

CHWP.TO 

HCVRTO 



0.002874 0.0101 

2.3323 

9.3599 



HEAT BALANCE: 





DHA DHH 

DHC 




4.263 0.079 

9.601 




COOLING AIR: 





75 T5 

T6 

HC 

ETA 


5.149 68.000 

69.63? 

10.244 

0.799 


PREFERRED VALUE OF DELT1 IS 15.62 UNDER 

THESE CONDITIONS 

SO THAT 72.62 = T2 = 

DPI + DELT1 

= o?*00 

+ 15.62. 


■a a a «• « «■ a j m p i j t D A T A ^ ^ * * * * * *' * 




CARD If PROCESS AIR 

I N : 




PCI P 01 PSA 

P1PSA 

T1 

DPI 

VI 

1.50 3*10 

14.70 

51.00 

46.00 

40*00 

CARD 2: H2 IN: 





V3SL T3 

DP3 

P3PSG 



1.000 70.000 

46 .000 

0.000 



CARD 3: MODULE: 





N I 

DELT1 

T5 

NFL Ac 


90. 4.88 

20.00 

40.00 

0 






INSUFFICIENT HYDROGEN. 4. 

33 SLPM 

IS REQUIRED 

* 

AVERAGE PROCESS AIR 

DEN POINT DEPRESSION OUT OF RANGE. 

ELECTROLYTE MOISTURE 

BALANCE IS 

MOT MAINTAINED. 



D-4 
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MODULE: 


T2 

66.000 

E 

0.359 

ESTACK 

32.338 

CURDEN 

20.000 

POWER 

157.808 

HEATLD 

1337.485 

OCON 

0.28901 

HC ON 
0.03641 

WPKOD 

0.32542 

TI 

1.3601 

TE 

49.4578 

C TRANS 
0.3931 

PROCESS AIR 
POl 
160.272 

IN: 

PI 

760.000 

P T W1 

7-929 

PN1 

590.299 

PW 1 
1.166 

FC1 

0.539 

RH1 

82.933 

VOl 

8.435 

VW1 

0.417 




PROCESS AIR 
P2 

760.00 

OUT: 

PC2 

0.41 

Pw'2 

10.15 

V2 

40.00 

VC2 

0.0214 

VW2 

0-5344 

VN2 

31.068 

V 02 
8.377 

PC 2 
0-146 

Yu 2 
1.493 

DP 2 
52.632 

RH2 

62.049 


H2 IN: 

V3 PW3 P 3 PH 3 VH3 I'H3 

3-531E-02 7.929 760.0 752.1 3.684K-04 3.494E-02 1.093E-02 


ANODE GAS OUT: 


DP 4 
49.316 

PW4 

8.981 

?C4 

0-393 

V H4 
-.081473 

VC4 

0.057563 

?4 

760.00 

Vis 4 

-.000286 - 

Fit 4 
.000799 

V4 

- . 0242 

V4SL 

-.6852' 

VH4SL 

-2.3073 

VC4SL 

1.0302 

FH4 

-.025485 

F4 

0.3668 

CHWRTO 

-15.4239 

KCVRTO 

-1.4154 



HEAT BALANCE 
DBA 
645.900 

DHH 

-.151 

DHC 

691.736 




COOLING AIR: 
v5 

62.883 

T5 

40.000 

T6 

50.219 

HC 

1.748 

ETA 

0.970 



PREFERRED VALUE OF DELil IS 23*63 UNDER THESE CONDITIONS 
SO THAT 69.65 = T2 = DPI 4 - DELT1 = 46.00 + 25. 63^ 


D-S 
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<hhhhmhhhh* I Npui DA TA*****^ # *** 
CARD 1* PROCESS AIR IN: 


PCI P01PSA 

P1PSA 

T1 

DPI 

VI 

10.00 3s 10 

14.70 

51.00 

46.00 

40.00 

CARD 2: H2 IN: 





V3SL T3 

DP3 

P3PSG 



20.000 70.000 

46.000 

0.000 



CARD 3: MODULE: 





N I 

DELT1 

T5 

NFL AG 


90. 4.88 

20.00 

90.00 

0 



INSUFFICIENTLY LOW COOLING AIR TEMPERATURE. 

AVERAGE PROCESS AIR DEW POINT DEPRESSION OUT OF RANGE. 
ELECTROLYTE MOISTURE BALANCE IS NOT MAINTAINED. 


MODULE: 


T2 E 

66.000 0.359 

ESTACK 

32.338 

CURDEN 

20.000 

POWER 
157. Q08 

HEATLD 

1337.485 

OC ON HCON 

0.28901 0.03641 

WPROD 

0.32542 

SI 

2.6488 

TE 

96.3200 

C TRANS 
0.7655 

PROCESS AIR IN: 





POl PI 

160.272 760.000 

pwi 

7.929 

PN 1 

581.799 

FW 1 

1-166 

FC1 

3.594 

RH 1 VOl 

82.933 3.435 

VW 1 

0.417 




PROCESS AIR OUT: 





P2 PC2 

760.00 7 088 

PW2 

10-14 

V2 

39.95 

VC2 

0.4142 

VW2 

0.5328 

VN2 ¥02 

30.621 8.377 

FC 2 

2.828 

FW2 

1.489 

DP2 

52.590 

KH 2 

61.953 

H2 IN: 





V3 PW3 

0.70G2 7.929 

P3 

760.0 

PH3 

752-1 

VW3 

7.368E-03 

VH3 

0.C988 

ANODE GAS OUT: 





DP 4 PW4 

49.295 8.957 

FC4 

0.766 

VH4 

0.582432 

VC4 

0.112106 

P4 

760.00 

VW4 FW4 

0.008233 0.023146 

V4 

0 .7028 

V4SL 

19.9039 

VH4SL 

16.4945 

VC4SL 

3.1748 

FH4 F4 

0.182185 0.9709 

CfisIRTO 

4.2019 

HCVRTO 

5.1954 




FH3 

0.2186 




Cife Systems. Jhc. 


HEM BALANCE: 


DHA 

645.900 

DHH 

“3.013 

DHC 

694.593 



COOLING AIRs 

V5 

T5 

T6 

HC 

ETA 

0 ° 0 0 0 

90.000 

0.000 

0.000 

0.000 


PREFERRED VALUE OF DELT1 IS 23.59 UNDER THESE CONDITIONS 
SO THAT 69.59 = T2 = DPI + DELT1 = 46.00 + 23.59. 

N j J U T DATA^*" JHt, ‘* HHKHHi ' 

CARD It PROCESS AIR IN: 


. PCI P01PSA 

PI PSA 

T1 

DPI 

VI 

3.00 3.10 

14.70 

51.00 

46.00 

20.00 

CARD 2; H2 IH: 





V3SL X3 

DP3 

P3PSG 



5.000 70.000 

46.000 

0.000 



CARD 3: MODULE; 





N I 

DELT1 

15 

NFL AG 


45. 9«76 

25.00 

51.00 

0 
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AVERAGE PROCESS AIR 

DEW POINT DEPRESSION 

OUT OF RANGE. 

RLE C TR OL Y TS M 0 1 STURE 

BALANCE IS 

NOT MAINTAINED. 


MODULE; 





T2 E 

ESTACK 

CURDaN 

POWER 

HEATLD 

71.000 0.271 

12.139 

40 *000 

118.960 

1470.307 

OC ON HCON 

ivPROD 

XI 

TS 

C TRANS 

0.28901 0.03641 

0.32542 

1.2349 

44.9043 

0.3569 

PROCESS AIR IN: 





POl PI 

Pwl 

PK1 

FBI 

FC1 

160.272 760.000 

7.929 

588.799 

0 .533 

0.539 

RH1 VOl 

VW1 




82.933 4.218 

0-209 




PROCESS AIR OUT: 





P2 PC2 

PW2 

V2 

VC2 

VS 2 

760.00 1-01 

12.36 

20.01 

0.0267 

0 .3255 

VM2 V02 

FC2 

FK2 

DP2 

RH2 

1 d . 49o 4*159 

0.182 

0.910 

58.057 

63 .'649 

H2 IN: 





V5 PW3 

P3 

PH3 

VS 3 

VH3 FH3 

0.1766 7.929 

760.0 

752.1 

1.842E-03 

0-1747 5 • 465 E® 02 


D-7 


£ife Systems, Juc. 


AHODE SAS OUT: 


DP4 

52.028 

PW4 

9.923 

FC4 

0.357 

VH4 

0.058296 

VC4 

0.052264 

P4 

760.00 

VW4 

0.001463 0 

FW4 

.004087 

V4 

0.1120 

V4SL 

3.1725 

VH4SL 

1.6509 

YC4SL 

1.4801 

FH4 

0*018235 

F4 

0.3792 

CHWRTO 

19.5713 

HCVRTO 

1.1154 



HEAT BALANCE 
DHA 
430.600 

0 

* 

DHH 

0.188 

DHC 

1039.519 




COOLING AIRS 
¥5 

299.394 

T5 

51.000 

T6 

54.225 

HC 

12.019 

ETA 

0.772 



PREFERRED VALUE OF DEL'i'l IS 29.06 UNDER THESE CONDITIONS 
SO THAT 75.06 = T2 = DPI + DELT1 = 46.00 + 29.06. 



